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Abstract. Internal diffusion-limited aggregation is a growth model 
based on random walk in TL'^ . We study how the shape of the aggre- 
gate depends on the law of the underlying walk, focusing on a family of 
walks in Tl? for which the limiting shape is a diamond. Certain of these 
walks — those with a directional bias toward the origin — have at most 
logarithmic fluctuations around the limiting shape. This contrasts with 
the simple random walk, where the limiting shape is a disk and the best 
known bound on the fluctuations, due to Lawler, is a power law. Our 
walks enjoy a uniform layering property which simplifles many of the 
proofs. 



1. Introduction and main results 

Internal diffusion-limited aggregation (internal DLA) is a growth model 
proposed by Diaconis and Fulton [DF91]. In the original model on Z"', parti- 
cles are released one by one from the origin o and perform simple symmetric 
discrete-time random walks. Starting from the set A(l) = {o}, the clusters 
A{i^V) for i > 1 are defined recursively by letting the i-th particle walk until 
it first visits a site not in A{%\ then adding this site to the cluster. Lawler, 
Bramson and Griffeath [LBG92] proved that in any dimension d > 2, the 
asymptotic shape of the cluster A(i) is a d- dimensional ball. Lawler [La95] 
subsequently showed that the fluctuations around a ball of radius r are at 
most of order r^/^ up to logarithmic corrections. Moore and Machta [MMOO] 
found experimentally that the fluctuations appear to be at most logarithmic 
in r, but there is still no rigorous bound to match their simulations. Other 
studies of internal DLA include [GQOO, BQR03, BB07, LP09b]. 

Here we investigate how the shape of an internal DLA cluster depends on 
the law of the underlying random walk. Perhaps surprisingly, small changes 
in the law can dramatically affect the limiting shape. Consider the walk 
in with the same law as simple random walk except on the x and y-axes, 
where steps toward the origin are reflected. For example, from a site (x,0) 
on the positive x-axis, the walk steps to 0) with probability 1/2 and to 

each of (x, ±1) with probability 1/4; see Figure 1. It follows from Theorem 1, 
below, that when we rescale the resulting internal DLA cluster A(i) to have 
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Figure 1. Example of a uniformly layered walk. The sites 
enclosed by the shaded area form the diamond V^. Only 
the transition probabilities from layer are shown. Open- 
headed arrows indicate transitions that take place with prob- 
ability 1/2; all the other transitions have probability 1/4. 

area 2, its asymptotic shape as i — oo is the diamond 

P = {(x,y) eR' : |x| + |y| < 1}. 

In fact, a rather large family of walks produce this diamond as their 
limiting shape. The key property shared by the walks we will consider is 
that their position at any time t is distributed as a mixture of uniform 
distributions on diamond layers. To define these walks, for A; > let 

jCk ■= {x^TI? : ||x|| = k} 

where for x = {xi,X2) we write ||x|| = \xi\ + \x2\- A uniformly layered 
walk is a discrete-time Markov chain on state space whose transition 
probabilities Q{x,y) satisfy 

(Ul) Q{x,y) =Oif\\y\\ > \\x\\ + l; 

(U2) For all A; > and all x G Ck, there exists y G JCk+i with Q{x, y) > 0; 
(U3) For all k,i>0 and all y,ze Ce, 

x^Ck xeCk 
In order to state our main results, let us now give a more precise descrip- 
tion of the aggregation rules. Set A{1) = {o}, and let Y^{t) (i = 1, 2, . . . ) be 
independent uniformly layered walks with the same law, started from the 
origin. For i > 1, define the stopping times cr* and the growing cluster A{i) 
recursively by setting 

a' = mm{t >0:Y'\t)^ A{i)} 

and 

A{i+1) = A{i)Li{Y\a')}. 
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Figure 2. Internal DLA clusters in based on the uni- 
formly layered walk with transition kernel pQin + ^Qout- 
Left: p = 0, walks are directed outward. Right: p = 1/2, 
walks have no directional bias. Each cluster is composed of 
■^350 = 24 5 701 particles. 

Now for any real number r > 0, let 

Vr := {x G : < r} . 

We call Dr the diamond of radius r in Z^. Note that Vr = -D[rj • For integer 
n > 0, we have V^ = Ufc=o Since = 4:k for A; > 1, the volume of Vn 
is 

Vn := #P„ = 2n(n + l) + l. 

Our first result says that the internal DLA cluster of Vn sites based on any 
uniformly layered walk is close to a diamond of radius n. 

Theorem 1. For any uniformly layered walk in Tl? , the internal DLA clus- 
ters A{vn) satisfy 

^ (^n-4V^n^ C A{vn) C ^?„+2o0n^ eventually^ = 1. 

Here and throughout this paper eventually means "for all but finitely 
many n." Likewise, we will write i.o. or infinitely often to abbreviate "for 
infinitely many n." 

Our proof of Theorem 1 in Section 5 follows the strategy of Lawler [La95] . 
The uniform layering property (U3) takes the place of the Green's function 
estimates used in that paper, and substantially simplifies some of the argu- 
ments. 

Within the family of uniformly layered walks, we study how the law of the 
walk affects the fluctuations of the internal DLA cluster around the limiting 
diamond shape. A natural walk to start with is the outward-directed layered 
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probabilities of the outward di- 
transition probabilities for the 
The origin o is near the lower- 



walk X{t) satisfying 

\\x(t + i)\\ = ||xm|| + i 



for all t. There is a unique such walk satisfying condition (U3) whose tran- 
sition probabilities are symmetric with respect to reflection about the axes. 
It is defined in the first quadrant by 

Qont(Xx,y),{x,y + 1)) = for = 1,2, . . . , (1.1) 

X + y + I 

Qout((x,y),(x + l,y)) = for x,y = 1,2,..., (1.2) 

X + y + 1 

and on the positive horizontal axis by 

Qout((x,0),(x,±l)) = for x = 1,2,..., (1.3) 

X + 1 

Qout((a:,0),(x + l,0)) = for x = 1,2,.... (1.4) 

X + \ 

In the other quadrants Qout is defined by reflection symmetry, and at the 
origin we set (5out(o, z) = 1/A for all 2: G with ||2;|| = 1. See Figure 3. 

Likewise one can construct a symmetric Markov kernel defining an inward 
directed random walk which remains uniformly distributed on diamond lay- 
ers. This kernel is defined in the first quadrant by 

Q\n{{x,y),{x,y -I)) = ^ for = 1,2,. . . , (1.5) 

' X + y — \ 

X — 1/2 

Qin{{x,y),{x -l,y)) = — for = 1,2,. . . , (1.6) 

X + y — 1 
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p = l/2 p = 'i/A 



Figure 4. Closeups of the boundary of the diamond. Fluc- 
tuations decrease as the directional bias of the walk tends 
from outward {p = 0) to inward {p = 1). 

and on the positive horizontal axis by 

Qin((2;,0),(x-1,0)) =1 for x = 1,2,.... (1.7) 

Again, the definition extends to the other quadrants by reflection symmetry, 
and is completed by making the origin an absorbing state: Qin(o, o) = 1. 
See Figure 3. 

We now choose a parameter p S [0, 1), let q = 1 — p and define the kernel 
Qp := pQin + qQout- The parameter p allows us to interpolate between a 
fully outward directed walk at p = and a fully inward directed walk at 
p = 1. 

Theorem 1 shows that the fluctuations around the limit shape are at 
most of order y/ n log n for the entire family of walks Qp. However, one 
may expect that the true size of the fluctuations depends on p. When p 
is large, particles tend to take a longer time to leave a diamond of given 
radius, affording them more opportunity to fill in unoccupied sites near the 
boundary of the cluster. Indeed, in simulations we find that the boundary 
becomes less ragged as p increases (Figure 4). Our next result shows that 
when p> 1/2, the boundary fluctuations are at most logarithmic in n. 

Theorem 2. For all p ^ {1/2,1), we have 

F (P„-6iog,.n C A{vn) C 2?„+6iog,.n eventually) = 1 

where the base of the logarithm is r = p/q. 
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Figure 5. Left: Example of a uniformly layered walk on 
the triangular lattice with hexagonal layers. Only transitions 
from a single (shaded) layer are shown. Open-headed arrows 
indicate transitions that take place with probability 1/2; all 
the other transitions have probability 1/4. Right: An inter- 
nal DLA cluster of 100 000 particles based on this uniformly 
layered walk. 

We believe that for all p G [0, 1/2) the boundary fluctuations are of or- 
der ^/n up to logarithmic corrections, and that therefore an abrupt change 
in the order of the fluctuations takes place at p = 1/2. At present, however, 
we are able to prove a lower bound on the order of fluctuations only in the 
case p = 0: 

Theorem 3. For p = we have 

^ K-(l^.)V 2(nloglog.)/3 ^ ^(^n) ^.0.) =1 V. > 

and 

Uniformly layered walks are closely related to the walks studied in [Du04, 
Ka07]. Indeed, the diamond shape of the layers does not play an important 
role in our arguments. A result similar to Theorem 1 will hold for walks 
satisfying (U1)-(U3) for other types of layers Ck, provided the cardinality 
#Ck grows at most polynomially in k. Figure 5 shows an example of a walk 
on the triangular lattice satisfying (U1)-(U3) for hexagonal layers. The 
resulting internal DLA clusters have the regular hexagon as their asymp- 
totic shape. Blachere and Brofferio [BB07] study internal DLA based on 
uniformly layered walks for which i^Ck grows exponentially, such as simple 
random walk on a regular tree. 

Given how sensitive the shape of an internal DLA cluster is to the law of 
the underlying walk, it is surprising how robust the shape is to other types of 
changes in the model. For example, the particles may perform deterministic 
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rotor-router walks instead of simple random walks. These walks depend 
on an initial choice of rotors at each site in Z'^, but for any such choice, 
the limiting shape is a ball. Another variant is the divisible sandpile model, 
which replaces the discrete particles by a continuous amount of mass at each 
lattice site. Its limiting shape is also a ball. These models are discussed 
in [LP09a]. 

The remainder of the paper is organized as follows. Section 2 explores 
the properties of uniformly layered walks, section 3 discusses an "abelian 
property" of internal DLA which is essential for the proof of Theorem 1, 
and section 4 collects the limit theorems we will use. Sections 5, 6 and 7 are 
devoted to the proofs of Theorems 1, 2 and 3, respectively. 



2. Uniformly layered walks 

Let {X(t)}t>o be a uniformly layered walk, that is, a walk on satisfying 
properties (U1)-(U3) of the introduction. Write i^k for the uniform measure 
on the sites of layer C^, and let denote the law of the walk started 
from X{0) ~ z^fc. Likewise, let F^; denote the law of the walk started from 
X{0) = X. Consider the stopping times 

:= mm{t > : X{t) = z] for z G Z^; 

Tfc := min{t > : X(t) G for /c > 0. 

The key to the diamond shape, as we shall see, is the fact that our random 
walks have the uniform distribution on diamond layers at all fixed times, 
and at the particular stopping times r^. The next lemma shows that un- 
der Ffc, conditionally on ||X(s)|| for s <t, the distribution oi X(t) is uniform 
on We remark that the fact that this conditional distribution de- 

pends only on and not on ||X(s)|| for s < t, implies that ||X(t)|| is 

a Markov chain under F/t; see [RP81]. 

Lemma 4. Fix k > 0. For all t > and all sequences of nonnegative 
integers k = ■ ■ ■ , i{t) satisfying £{s -|- 1) < £{s) + 1 for s = 0, . . . , t — 1, 
we have for all z G 

Fk{X{t) = z I \\X{s)\\ = £{s), 0<s<t) = 

= Vk{X{t) = z\ \\X{t)\\=£{t)). 

Proof. We prove the first equality by induction on t. The base case t = is 
immediate. Write 



£t = {\\X{s)\\=£{s),Q<s<t]. 
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By the Markov property and the inductive hypothesis, we have for t > 1 
and any y G C^^t) 



Fk{X{t) = y, 8t) = ^k{X{t) = y, X{t -l) = x, £t- 

= Yl Q{x,y)-'Pk{x{t-i) =x,£t^i 



e(t-i) 



Yl Q{x,y)-^^J—-lPk{£t-i). 



By property (U3), the right side does not depend on the choice of y G J~-i{t)- 
It follows that 



By induction this holds for all t > and all sequences ^(0), . . . , £{t). There- 
fore, for fixed £{t) and z G 

Ffe {X{t) = z)= Y (^(^) = ^' II = ^(^) ^ 

£(o),...,e{t-i) 

= #r' E Pk{\\x{s)\\=iis)ys<t) 

^ i{o),...,e{t-i) 



^ P,.(||X(t)||=^(t)) 



which implies 



Fk{x{t) = z I \\x{t)\\ = m) = ^7^- □ 



As a consequence of Lemma 4, our random walks have the uniform dis- 
tribution on layer i at the stopping time T£. 

Lemma 5. Fix integers < k < i. Then 

^k{X{Ti) = z) = = ^ for every z G Ci>. 

Proof. Note that property (U2) and Lemma 4 imply < cxd almost surely. 
For t > we have 

{n = t]= (J {||x(s)|| =4,o<s<t}, 

£o,-/t 

where the union is over all sequences of nonnegative integers ^Oj^i; • • • i^t 
with £q = k and £t = such that ig+i < + 1 and is i for all s = 
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0, 1, . . . , t — 1. Writing £eo,...,it foi' the disjoint events in this union, it follows 
that 

Fk{X{Te) = z) = Y^ Ffc(X(t) =z,Te = t) 

t>0 lo,.../t 

= E E ^k{X{t) = z\ k{£i,,...,t,). 

t>0 to,.../t 

Since YltX)^ k{Tii = = 1, the result follows from Lemma 4. □ 

The previous lemmas show that one can view our random walks as walks 
that move from layer to layer on the lattice, while remaining uniformly dis- 
tributed on these layers. This idea can be formalized in terms of an inter- 
twining relation between our two-dimensional walks and a one-dimensional 
walk that describes the transitions between layers, an idea explored in [Du04, 
Ka07] for closely related random walks in wedges. This approach is partic- 
ularly useful for computing properties of the Green's function. 

Next we calculate some hitting probabilities for the walk with transition 
kernel Qp = p Qin + Q Qcmt defined in the introduction; we will use these 
in the proof of Theorem 2. We start with the probability of visiting the 
origin before leaving the diamond of radius n. By the definition of Qp, this 
probability depends only on the layer on which the walk is started, not on the 
particular starting point on that layer. That is, if < £ < n, then Fa;(To < 
Tn) = F^(to < T„) for all X S Ce, since at every site except the origin, 
the probability to move inward is p and the probability to move outward 
is q. This leads to the following well-known gambler's ruin calculation (see, 
e.g., [Bi!).>, §7]). 

Lemma 6. Let < i < n and x £ Cg. If p ^ q, then 

r"- — 1 

where r = p/q. If p = q = 1/2, then 

n — i 

Fx(To < Tn) = Fi{To < Tn) = . 

n 

Next we bound the probability that the inward- biased walk (p > 1/2) 
exits the diamond T>n-i before hitting a given site z £ T>n-i- 

Lemma 7. Write r = p/q. For p G (1/2, 1), if z £ for < k < n, then 

Fo(r. > Tn) < {Ak - l)r'^-". 

Proof. Let Tq = and for i > 1 consider the stopping times 

Ui = mm{t > : X{t) G 

Ti = mm{t > Ui : X{t) = o}. 
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Let M = max{z : Ui < Tn}- For any integer m>l and any xi, . . . , Xm £ Ck-, 
we have by the strong Markov property 

Fo(M = m, X{Ul) = XI, . . . , X{Um) = Xm) 

m— 1 

= J] [Fo(X(rfc) = Xi) ¥^^{to < Tn)] ■ Fo(X(rfc) = x^) F^,„(r„ < r^). 

i=l 

By Lemma 5, Fo(X(Tfc) = Xj) = 1/4A; for each Xj G £fc. Moreover, by 
Lemma 6 we have for any x S £fc 

F.(t„ < to) = 4^ < 
r" — 1 

where we have used the fact that r = p/q > 1. Hence 

Fo(M = m, XiUi) i^z'ii<m) < r'^"" f 1 - ^ 

Since the event {r^ > r„} is contained in the event {X{Ui) ^ z'ii < M}, 
we conclude that 



Fo(r, > r„) = ^ Fo(M = m, > t„) 

m>l 

< ^ Fo(M = m, X(;7i) / z Vi < m) 

m>l 

<E'-'-"(i-s)" 

m>l ^ ^ 

= (4A; - l)r'=-". □ 
3. AbELIAN PROPERTY 

In this section we discuss an important property of internal DLA discov- 
ered by Diaconis and Fulton [DF91, Theorem 4.1], which gives some freedom 
in how the clusters A{i) are constructed. We will use this property in the 
proof of Theorem 1. It was also used in [La95]. Instead of performing i 
random walks one at a time in sequence, start with i particles at the origin. 
At each time step, choose a site occupied by more than one particle, and 
let one particle take a single random walk step from that site. The abelian 
property says that regardless of these choices, the final set of i occupied sites 
has the same distribution as the cluster A{i). 

This property is not dependent on the law of the random walk, and in fact 
holds deterministically in a certain sense. Suppose that at each site x G 
we place an infinite stack of cards, each labeled by a site in Z^. A legal 
move consists of choosing a site x which has at least two particles, burning 
the top card at x, and then moving one particle from x to the site labeled 
by the card just burned. A finite sequence of legal moves is complete if it 
results in a configuration in which each site has at most one particle. 
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Lemma 8 (Abelian property). For any initial configuration of particles 
on TI? , if there is a complete finite sequence of legal moves, then any sequence 
of legal moves is finite, and any complete sequence of legal moves results in 
the same final configuration. 

In our setting, the cards in the stack at x have i.i.d. labels with distri- 
bution Q{x, ■). Starting with i particles at the origin, one way to construct 
a complete sequence of legal moves is to let each particle in turn perform 
a random walk until reaching an unoccupied site. The resulting set of oc- 
cupied sites is the internal DLA cluster A{i). By the abelian property, any 
other complete sequence of legal moves yields the same cluster A{i). 

For the proof of Theorem 1 , it will be useful to define generalized internal 
DLA clusters for which not all walks start at the origin. Given a (possi- 
bly random) sequence xi,X2, . . . G Z^, we define the clusters A{xi, . . . ,Xi) 
recursively by setting A{xi) = {xi}, and 

A{xi, Xi+i) = A{xi, ...,Xi)U {Y\a')}, i > 1, 

where the are independent uniformly layered walks started from Y^(0) = 
Xj+i, and 

a' = mm{t > : Y\t) ^ A{xi, Xi)}. 

When xi = ■ ■ ■ = Xi = o we recover the usual cluster A(i). 

The next lemma gives conditions under which two such generalized clus- 
ters can be coupled so that one is contained in the other. Let xi, . . . ,Xr and 
yi, . . . ,ys be random points in Z^. For z ^TI? ^ let 

A^. = #{« < r : = z] 
Nz = #{j <s : yj = z} 

and consider the event 

Lemma 9 (Monotonicity) . There exists a random set A' with the same 
distribution as A{yi, . . . , y^), such that £ C {A{xi, . . . , Xr) C A'^. 

The proof follows directly from the abelian property: since the distribu- 
tion of A{yi, . . . ,ys) does not depend on the ordering of the points j/i, . . . , y^, 
we can take 

U(y'i,---,ys) on£ 

\A{yi,...,ys) on£'' 

where y[, . . . ,y'g is a (random) permutation of yi , . . . , ?/s such that y'^ = xi 
for all i < r. 
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4. Sums of independent random variables 

We collect here a few standard results about sums of independent random 
variables. First we consider large deviation bounds for sums of independent 
indicators, which we will use several times in the proofs of Theorems 1 and 2. 
Let S he a finite sum of independent indicator random variables. We start 
with simple Chernoff-type bounds based on the inequality 

F(5 > 6) < e-^^'E (e*-^) . 

There are various ways to give an upper bound on the right side when the 
summands of S are i.i.d. indicators; see for example [AS92, Appendix A]. 
These bounds extend to the case of independent but not necessarily identi- 
cally distributed indicators by an application of Jensen's inequality, leading 
to the following bounds " ', Theorems 1 and 2]: 

Lemma 10 (Chernoff bounds). Let S be a finite sum of independent indi- 
cator random variables. For all b >0, 

. 1 b^ 

F{S > + 6) < exp 



F{S < - 6) < exp 



2E5 + 5/3 
2ES' 



Next we consider limit theorems for sums of independent random vari- 
ables, which we will use in the proof of Theorem 3. For {X„}„>i a sequence 
of independent random variables satisfying E|Xj|^ < cxd, we define 



Bn = Yl Var(X,), (4.1) 

l<j<n 

L„ = 5-3/2 ^ E|Xi-EXi|3. (4.2) 

i<j<?i 

It is well known that the partial sums 

Sn= (4-3) 

l<i<n 

satisfy the Central Limit Theorem when Ln — > 0; this is a special case of 
Lyapunov's condition. We are interested in the rate of convergence. Let 

F (Sn - E5„ < x^) - $(x) , (4.4) 



A„ = sup 

where ^ is the standard normal distribution function. Esseen's inequality 
(see [Es45, Introduction, equation (6)] and [PSOO, Chapter I]) gives a bound 
on A„ in terms of L„. This bound can be used to verify the conditions given 
by Petrov [PeOfi, Theorem 1] (see also [PSOCI, Chapter I]), under which the 
partial sums Sn satisfy the Law of the Iterated Logarithm. 
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Lemma 11 (Esseen's inequality). Let Xi, . . . ,Xn be independent and such 
that < oo, and define Bn, Ln, Sn and hy (4.1)-(4.4). Then 

A„ < 7.5 • Ln. 

Lemma 12 (Petrov's theorem). Let {Xj}j>i he a sequence of independent 
random variables with finite variances, and define Bn, Sn and A„ hy (4.1), 
(4.3) and (4.4). //, as n ^ oo, 

-B„ ^ oo, ^ 1 0.'^'^ ^"^^(n R w+s ] for some 5 > f), 

then 

TO (v ~ ^^"^ 1^ 1 

F I iim sup —^^^^=^^^= = 1 = 1 . 

V n^oo ^/2Bnloglog Bn ) 



5. Proof of Theorem 1 

We control the growth of the cluster A{i) by relating it to two modified 
growth processes, the stopped process S{i) and the extended process E{i). In 
the stopped process, particles stop walking when they hit layer even if 
they have not yet found an unoccupied site. More formally, let S'(l) = {o}, 
and define the stopping times and clusters S{i) for i > 1 recursively by 

4 = min{t > : Y\t) G £„ U 5(i)'} 

and 

5(i + l) = 5(i)U{r(4)}. (5.1) 

Here Y^it) for i = 1,2, . . . are independent uniformly layered walks started 
from the origin in Tl? ^ all having the same law. Note that S{i + 1) = S{i) 
on the event that the walk hits layer Ln before exiting the cluster S{i). 
By the abelian property. Lemma 8, we have S{i) C A{i). Indeed, Aii) can 
be obtained from S{i) by letting all but one of the particles stopped at each 
site in continue walking until reaching an unoccupied site. 

The extended process £'(«) is defined by starting with every site in the 
diamond occupied, and letting each of i additional particles in turn 
walk from the origin until reaching an unoccupied site. More formally, let 
-E(O) = T>ni and for z > define 

4 = min{t > : y^"+Xt) ^ E{i)] 

and 

E(f + i) = E(f)u{y^"+^(4)}. 

An outline of the proof of Theorem 1 runs as follows. We first show in 
Lemma 13 that the stopped cluster S{vn) contains a large diamond with 
high probability. Since the stopped cluster is contained in A{vn), the inner 
bound of Theorem 1 follows. The proof of the outer bound proceeds in three 
steps. Lemma 14 shows that the particles that stop in layer during the 
stopped process cannot be too bunched up at any single site z G £„. We 
then use this to aiiguc in LGinma 15 that with, high probabihty, J\.(^Vfij is 
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contained in a suitable extended cluster E{m). Finally, Lemma 16 shows 
that this extended cluster is contained in a slightly larger diamond. 

A notable feature of the argument (also present in [La95]) is that the 
proof of the outer bound relies on the inner bound: Lemma 13 is used in the 
proof of Lemma 14. This dependence is responsible for the larger constant 
in the outer bound of Theorem 1. It would be interesting to know whether 
this asymmetry is merely an artifact of the proof, or whether the outer 
fluctuations are really larger than the inner fluctuations. 

We introduce an auxiliary collection of walks that will appear in the 
proofs. Let {Y^{t) : x £ Z^} be independent walks with the same transi- 
tion probabilities as Y^, which are independent of the Y^, and which start 
from y^(0) = X. Now for i = 1, — 1 define 



Y'{t) for < i < cr^, 

y^'('^s)(t-cjjj) fort>4. 



Note that replacing the walks Y'^ with X* in (5.1) has no effect on the 
clusters S{i). Finally, for i > Vn set X*(t) = Y^{t) for all t >0. 

We associate the following stopping times with the auxiliary walks Y^(t): 

rf := mm{t > : Y'^it) = z} for z G Z^; 

:= mm{t > : y^(t) G Ck} for A; > 0. 

Likewise, let 

ri := mm{t > : X\t) = z) for z £ Z^; 

4 := mm{t > : X\t) £ Ck} for k>0. 

Lemma 13. There exists uq such that for all uniformly layered walks and 
all n > no 



^ (^n-4vCT^ ^ SM) < 6n-2. (5.2) 



Remark. To avoid referring to too many unimportant constants, for the rest 
of this section we will take the phrase "for sufficiently large n," and its 
variants, to mean that a single bound on n applies to all uniformly layered 
walks. 

Proof. For z G I'n-i, write 

v„-l 
i=l 

for the event that the site z does not belong to the stopped cluster S{vn)- 
We want to show that P(^Sz{vn)) is very small when z is taken too deep 
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inside Vn. To this end, let £ = ||z||, and consider the random variables 











^ ^0<l<Vn 




Lz 


^ ^0<l<Vn 




Then £z{vn) = {N^ = 


0}. Since > - 


Lz, we have for 


number a 








P{Nz = 0) < F(M^ < a 


or Lz > a) 




< F(M^ < a 


) + ¥{Lz>a). 



(5.3) 

Our choice of a will be made below. Note that Mz is a sum of i.i.d. indicator 
random variables, and by Lemma 5, 

EM, = 2n(n + 1) Fo(X(t,) =z) = \ (5.4) 

The summands of Lz are not independent. Following [LBG92], however, we 
can dominate Lz by a sum of independent indicators as follows. By prop- 
erty (Ul), a uniformly layered walk cannot exit the diamond P^-i without 
passing through layer ££, so the event {ct^ < t| = t}} is contained in the 
event {X*(4) & T^i-i}- Hence 

Lz= Y: l{x>j,)eP,_„rf>^)=rf(-^)} 

0<j<Dn 

where we have used the fact that the locations X*(cj^) inside D^-i where 
particles attach to the cluster are distinct. Note that Lz is a sum of inde- 
pendent indicator random variables. To compute its expectation, note that 
for every < A; < ^, by Lemma 5 

^x{X{n) = z)= AkFk{X{Te) = z) = j, 

hence 

^Lz = Yj = '-^- (5-5) 

k=l 

Now set a = ^(EM, + EL,), and let 



EM, - EL, n 



2 ff2 



b= ^ ^ " > 



4£ 
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where the inequahty follows from (5.4) and (5.5). Since a = lEM^ — b 
+ b, we have by Lemma 10 

/I 62 \ / 1 62 

^{Lz >a)< exp -- — . < exp 



2 EL, + 6/3/ " " V 2 EM; 



and 



F{Mz <a)< exp 
< exp 



1 62 



2 EMz 

1 (^2 - £2)2 2£ 



2 16^2 n{n + 1) 

< exp 



1 (n2 



16 

where in the last line we have used i < n — 1. Since < L,, we obtain 
from (5.3) 

F{£z{vn)) < F(M, < a) + F{Lz > a) 

1 (n2 ^2^2 

< 2 exp ' 



16 

Writing £ = n — p, with p > |'4-y/nTogn] , we obtain for sufficiently large n 

1 p2(2n-p)2^ 



F{£z{vn)) < 2 exp 



16 



< 2 exp ( + y 



< 3n-^ 

We conclude that for n sufficiently large 

^^■'^n — 4Vi log n 

Turning to the outer bound of Theorem 1, the first step is to bound the 
number 

A^, := y 1{4 = tI] (5.6) 

of particles stopping at each site z G £„ in the course of the stopped process. 
To get a rough idea of the order of Nz, note that according to Lemma 13, 
with high probability, at least t^„_4^ n iog„ of the u„ particles find an occupied 
site before hitting layer £„. The number of particles remaining is of order 
n^/2^1og n. If these remaining particles were spread evenly over then 
there would be order yjn \ogn particles at each site z S £„. The following 
lemma shows that with high probability, all of the Nz are at most of this 
order. 
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Lemma 14. // n is sufficiently large, then 

^( U {^^ > 7x/nlogn} j < 13n-^/^. 
Proof. For z ^ Cn, define 

SO that Nz = Mz — Lz ■ Write rj = \/ n log n and p = [4//] , and let 

Lz= Yl li^' = ^n}- 

yeVn-p-lo} 

Note that Lz < Lz on the event {D^-p C Therefore, 

< Y nMz -Lz>7v) + F [V^.,^^ € S{vn) 

To obtain a bound on V[Mz — Lz > Jrj), note that 

n + 1 

EM, = 2n(n + 1) Fo(X(t„) = z) = ^— . 
Moreover, by Lemma 5 

k 
n 



Yp{ry = Ty) = 4kFk{Tz = Tn) 



hence 

n—p 



~ n + 1 , 1) 

EL, = > - = — p+ — . 

^-^ n 2 In 

k=l 

In particular, EM, — EL, < p — \ < Arj for large enough n, so that 

F(M, - L, > 7r/) < F(M, - L, > EM, - EL, + 3t/) 

< F(M, > EM, + §r/ or L, < EL, - §7?) 

< F(M, > EM, + |r/) +F(L, < EL, - |r/). 

By Lemma 10, 

P(4 < EL, - i,) < exp (4^) < exp = 
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Likewise, for sufficiently large n 

F{M^ > EM^ + Irj) < exp 



exp 



1 (3r?/2)" 
2EM2 + r?/2 
9 n log n 



4 n + 1 + y/ n log n 



Combining (5.7), (5.8) and Lemma 13 yields for sufficiently large n 



□ 



Given random sets A^B C Z^, we write A = i? to mean that A and i? 
have the same distribution. 



Lemma 15. Let m = |"29n-v/nTogn] . For all sufficiently large n, there exist 

random sets A' = A{vn) and E' = E{m) such that 

F(^' ^ E') < Un~^/^. 

Proof. By the abelian property. Lemma 8, we can obtain A{vn) from the 
stopped cluster S{vn) by starting particles at each z £ Cn, and letting 
all but one of them walk until finding an unoccupied site. More formally, 
let xi = and Xj+i = Y^{a^g) for < i < Vn- Then 



#{i <Vn : Xi = z} = < 



Nz, z e Cn 

1, Z £ S{Vn) - Cn 

0, else 



and 



A{Vn) = A{xi, . . .,XyJ. 

To build up the extended cluster E{m) in a similar fashion, let s = Vn+m, 
and let yi, . . . , i/s G be such that {yi, . . . , y^,,^} = Vn, and 

= y-"+-i«"+-i), i = 1,2, . . . ,m. 

By Lemma 8, we have 

E{m) = A{yi, . . .,ys). 

For each z £ Cn, let 

Nz = y 1 {tJ"+^ = r^"+n 

be the number of extended particles that first hit layer Cn at z. Then 





Nz, 


z e Cn 


#{i < s : yi = z} = < 


1, 


Z£Vn 




.0' 


else. 
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Now let A' = A{xi, . . . , Xv„) and consider the event 

f = fl {N,<N,}. 

zeC„ 

By Lemma 9, on the event £ there exists a random set E' = A{yi, . . . , t/s) 
such that A' C E'. Therefore, to finish the proof it suffices to show that 
W(£^) < 14n~^/^. Note that Nz is a sum of independent indicators, and 

^ m 29 



where rj := ^/rlTogn. Setting 6 = in Lemma fO yields for sufficiently 
large n 

F < 7'n] < exp f-^ = exp ( — \/n log ri] < ^n"^/'*, 

hence by Lemma f4 

F{£^)<f([J{Nz>77j or iV^ < Tr?} j 

< IJ {iV, > 7r?}^ + J] F(iV, < 77]) 

< Un-^/^. □ 

To finish the argument it remains to show that with high probability, the 
extended cluster E{m) is contained in a slightly larger diamond. Here we 
follow the strategy used in the proof of the outer bound in [LBG92]. 

Lemma 16. Let m = |"29n-v/nTogn] . For all sufficiently large n, 

F [E{m) t P„+20v/^) < 
Proof. For j,k > 1, let 

ZkU) = #{E{j) n Cn+k) 

and let ^kU) = lEZfc(j). Then fikU) the expected number of particles 
that have attached to the cluster in layer Cn+k after the first j extended 
particles have aggregated. Note that 

l^kii + 1) - l^kii) = F(y^"+^+Hai+i) G Cn+k). 

By property (Ul), in order for the (i + 1)*^ extended particle to attach to 
the cluster in layer Cn+k^ it must be inside the cluster E{i) when it first 
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reaches layer Cn+k-i- Therefore, by Lemma 5, 

+ 1) - ^.(i) < F(i^^"+*+'(CtV) e m) 



1 /.^ ^ ^fc-l(0 



4{n + k-l) ^^-'^ ' - An 
Since /ifc(O) = 0, summing over i yields 

1 



An 

i=l 



Since < j and ^^^^ ^ ^ j'^/k, we obtain by induction on /c 



4n J kl \4nk 

where in the last equality we have used the fact that kl > k^e~^. Since 
29e/80 < 1, setting j = m and k = [20-v/n log nj we obtain 



P -, \ fc+i 

I 29n V w log n 



4n • 20\/n log n / 

for sufficiently large n. To complete the proof, note that 

F{E{m) Vn+k) = F(^fc+i(m) > 1) < /ife+i(m). □ 

Proof of Theorem 1. Write r/ = y^nlogn. Since S{vn) C A{vn)^ we have by 
Lemma 13 

SiVn)) < OO. 

n>l n>l 

Likewise, by Lemmas 15 and 16 

Y,^{A{Vn) t Vn+2Qri) <Y,^{^M t S(m)) + ^ F (^(m) (/L Pn+20r,) 
n>l n>l n>l 

< OO. 

By Borel-Cantelli we obtain Theorem 1. □ 

6. The inward directed case 
Proof of Theorem 2. Write i = n — [6 log^ n] , and denote by 

0<i<Vn Z&Vg 

the event that each of the first f „ — 1 walks visits every site z £ before 
hitting layer Cn- Since < f„, at least one of the first f n — 1 particles 

must exit Vn-i before aggregating to the cluster: a' > for some i < Vn- 
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On the event An, this particle visits every site z G Vi before aggregating to 
the cluster, so Vi C A{i) C A{vn)- Hence 

0<i<Dn Z&Ve 

By Lemma 7, 

i 

F{Ve (f_ A{vn)) < 2n{n + 1) ^4A;(4fc - l)r''-'^ 

fc=i 

ij^^ I 1 

< 32n^(n + 1)— — 

< ^^n^(n + 1) • n~^, 
r — 1 

and by Borel-Cantelli we conclude that F(I'£ C A{vn) eventually) = 1. 
Likewise, writing m = n + [Olog^. nj , let 

o<i<vn zev 

be the event that each of the first u„ — 1 walks visits every site z G T>n 
before hitting layer Cm- Since the occupied cluster A{vn — 1) has cardinality 
f „ — 1 = #^?n — 1, there is at least one site z £ P„ belonging to A{vn — l^- 
On the event each of the first v„ — 1 particles visits z before hitting 
layer Cm, so 

(T"<Ti<Tl, i = l,...,Vn-l. 

Therefore, 

0<i<'!;„ z€Vn 

By Lemma 7, 



j.k-m 



F {A{vn) <t- T^m) < 2n{n + 1) X] ^^(4^ - 1) 

k=l 

< 32n^(n + l) — , — 

< -^^^n^(n + 1) • n"*^, 
r — 1 

and by Borel-Cantelli we conclude that W^A^Vn) C "Dm eventually) = 1. □ 
7. The outward directed case 

To prove Theorem 3 we make use of a specific property of the uniformly 
layered walks for p = 0. Recall that these walks have transition kernel Qout- 
By (1.1)-(1.4), such a walk can only reach the site (m, 0) for m > 1 by 
visiting the sites (0, 0), (1, 0), ... , (m, 0) in turn. We can use this to find the 
exact growth rate of the clusters A{i) along the x-axis. 
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Suppose that we count time according to the number of particles we have 
added to the growing cluster, and for m > 1 set 

Tjn := min{n > : (m, 0) € A{n + 1)}. 

Then we can interpret Tm as the time it takes before the site (m, 0) becomes 
occupied. The following lemma gives the exact order of the fluctuations 
in Tm as m — > oo. 

Lemma 17. For p = we have that 

Tm — 2m{m + 1) \ 



F limsup 



and 



m^oo 32(m^ log log m) /3 / 

r„-2,.(.,+ i) \ , 

ym^oc Y/32(m3loglogm)/3 J 

Proof. Set Xi = Ti and Xm = Tm — Tm-i for m > 1. Consider the aggre- 
gate at time Tm-i when (m — 1,0) gets occupied. Since a walk must follow 
the X-axis to reach the site (m — 1,0), we know that at time Tm~i all sites 
{(z,0) : i = 0, 1, . . . , m — 1} are occupied and all sites {(^,0) : i > m} are 
vacant. Now consider the additional time Xm = Tm — Tm-i taken before the 
site (m, 0) becomes occupied. Each walk visits (m, 0) if and only if it passes 
through the sites (1,0), (2,0),..., (m, 0) during the first m steps, which hap- 
pens with probability l/4m. Thus Xm has the geometric distribution with 
parameter l/4m. Moreover, the Xi are independent. Hence Tm is a sum of 
independent geometric random variables Xi. 

Since EXj = 4i, Var = IGz^ - Ai and EXf = 384%^ - 96i^ + 4i, 

16 

Bm= Y] VarXi = —m^ + 0{m^) 

o 

l<i<m 

and 

E(|Xi-EXi|3) < (EX3 + (EXi)3) =0(m4). 

l<i<m l<i<m 

By Lemma 11, Am = 0{m~^^'^), which shows that Petrov's conditions of 
Lemma 12 are satisfied. Therefore, 

(Tm — ^Tm \ 
lim sup - = 1 = 1 . 

■m^oo V2-Bm log log / 

Since ET^ = X^^ILi 4i = 2m(m-|-l) and Bm = 16m^/3 -|- O(m^), this proves 
the first statement in Lemma 17. The second statement is obtained by 
applying Lemma 12 to —Tm = J2^i{—Xi). □ 
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Proof of Theorem 3. Fix e > 0, set r/ := y^2(nloglogn)/3 and let p = 
[(1 — e)r/] . If we write m = n — p, then 

2m{m + 1) = 2n(n + 1) — 4np + o(n), 

■\/ 32 (m^ log log m) /3 = 4?7,r? + o(n^/^ log log n) . 

Hence, setting m = n — p in Lemma 17 gives 

pflimsup^"---^"'"+" + ^"^ = lUL 

V n^oo 4nr/ y 

Since {(n — 0) = {T„_p > f„, — 1} and Vn — 1 = 2n(n + 1), this 

implies 

F{{n-p,0)^A{vn) i.o.) = 1. 
Likewise, setting m = n + p in Lemma 17 gives 



p(limmf^"^'--'^"'"+"-''"^ = -lUl. 



hence 



F{{n + p,0) e A{vn) i.o.) = 1. □ 

8. Concluding Remarks 

Lawler, Bramson and GrifFeath [LBG92] discovered a key property of the 
Euclidean ball that characterizes it as the limiting shape of internal DLA 
clusters based on simple random walk in Z'^: for simple random walk killed 
on exiting the ball, any point z sufficiently far from the boundary of the ball 
is visited more often in expectation by a walk started at the origin than by 
a walk started at a uniform point in the ball. Uniformly layered walks have 
an analogous property with respect to the diamond: the Green's function 
g(y, •) for a walk started at y and killed on exiting satisfies 

g{o,z) > g{y,z) 

for all z S T>n. Indeed, both the walk started at o and the walk started at a 
uniform point in T>n are uniformly distributed on layer at the time t\\z\\ 
when they first hit this layer, so the expected number of visits to z after 
time T||2|j is the same for both walks. The inequality comes from the fact 
that a walk started at the origin must hit layer before exiting P„. 

We conclude with two questions. The first concerns uniformly layered 
walks started from a point other than the origin. Figure 6 shows internal 
DLA clusters for six different starting points in the first quadrant of Z^. 
These clusters are all contained in the first quadrant. Our simulations in- 
dicate that a limiting shape exists for each starting point, and that no two 
starting points have the same limiting shape; but we do not know of any 
explicit characterization of the shapes arising in this way. 

The second question is, do there exist walks with bounded increments 
having uniform harmonic measure on spheres in Z*^ for d > 3? 
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(2,2) (3,2) 




(1,1) (2,1) (3,1) 

Figure 6. Internal DLA clusters in the first quadrant of H? 
based on the outward-directed layered walk Qout started from 
a point other than the origin. For example, the cluster on 
the lower left is formed from 405 900 particles started at the 
point (1, 1). 
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